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Abstract 

Big-bang nucleosynthesis (BBN) describes the production of the lightest nuclides via a dynamic 
interplay among the four fundamental forces during the first seconds of cosmic time. We briefly 
overview the essentials of this physics, and present new calculations of light element abundances 
through 6 Li and 'Li, with updated nuclear reactions and uncertainties including those in the 
neutron lifetime. We provide fits to these results as a function of baryon density and of the number 
of neutrino flavors, N u . We review recent developments in BBN, particularly new, precision Planck 
cosmic microwave background (CMB) measurements that now probe the baryon density, helium 
content, and the effective number of degrees of freedom, N e g. These measurements allow for a 
tight test of BBN and of cosmology using CMB data alone. Our likelihood analysis convolves 
the 2015 Planck data chains with our BBN output and observational data. Adding astronomical 
measurements of light elements strengthens the power of BBN. We include a new determination 
of the primordial helium abundance in our likelihood analysis. New D/H observations are now 
more precise than the corresponding theoretical predictions, and are consistent with the Standard 
Model and the Planck baryon density. Moreover, D/H now provides a tight measurement of N v 
when combined with the CMB baryon density, and provides a 2cr upper limit N v < 3.2. The 
new precision of the CMB and of D/H observations together leave D/H predictions as the largest 
source of uncertainties. Future improvement in BBN calculations will therefore rely on improved 
nuclear cross section data. In contrast with D/H and 4 He, 7 Li predictions continue to disagree 
with observations, perhaps pointing to new physics. We conclude with a look at future directions 
including key nuclear reactions, astronomical observations, and theoretical issues. 
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I. INTRODUCTION 


Big bang nucleosynthesis (BBN) is one of few probes of the very early Universe with 
direct experimental or observational consequences In the context of the Standard 

Models of cosmology and of nuclear and particle physics, BBN is an effectively parameter- 
free theory fSJ. Namely, standard BBN (SBBN) assumes spacetime characterized by General 
Relativity and the ACDM cosmology, and microphysics characterized by Standard Model 
particle content and interactions, with three light neutrino species, with negligible effects 
due to dark matter and dark energy. 

In SBBN, the abundances of the four light nuclei are usually parameterized by the baryon- 
to-photon ratio = nb/n 7 , or equivalently the present baryon density, i\K 2 = c < 7 ,. This quan¬ 
tity has been fixed by a series of precise measurements of microwave background anisotropies, 
most recently by Planck yielding r] = 6.10 ± 0.04 [5]. Thus the success or failure of SBBN 
rests solely on the comparison of theoretical predictions with observational determinations. 

While precise predictions from SBBN are feasible, they rely on well-measured cross sec¬ 
tions and a well-measured neutron lifetime. Indeed, even prior to the WMAP e ra, theoretical 
predictions for D, 3 He, and 4 He were reasonably accurate, however, uncertainties in nuclear 
cross sections leading to 'Be and 'Li were relatively large. Many modern analyses of nuclear 
rates for BBN were based on the NACRE compilation [7] and recent BBN calculations by 
several groups are in good agreement MZJ. Recent remeasurements of the 3 He(a, 7 ) 7 Be 
cross section p 2 j did improve the theoretical accuracy of the prediction but exacerbated the 
discrepancy between theory and observation [19j. Very recently, the NACRE collaboration 
has issued an update (NACRE-II) of its nuclear rate tabulation [20]. These were first used 
in [ 2 Tj] and we incorporate the new rates in the results discussed below. 

The neutron mean life has had a rather sordid history. Because it scales the weak inter¬ 
action rates between n p, the neutron mean life controls the neutron-to-proton ratio at 
freeze-out and directly affects the 4 He abundance (and the other light elements to a lesser 
extent). The value r n = 918 ± 14 s reported by Cristensen et ah in 1972 j 22 | dominated 
the weighted mean for the accepted value through mid 1980s. Despite the low value of 877 
± 8 s reported by Bondarenko et al. in 1978 [23], the ‘accepted’ mean value (as reported 
in the “Review of Particle Physics”) remained high, and the high value was reinforced by 
a measurement by Byrne et ah [24] hr 1980 of 937 ± 18 s. The range 877 s - 937 s was 
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used by Olive et al. [25] to explore the sensitivity of BBN predictions to this apparently 
uncertain quantity treated then as an uncertain input parameter to BBN calculations (along 
with the number of light neutrino flavors N v and the baryon-to-photon ration, rf). In the 
late 1980’s a number of lower measurements began to surface, and in 1989 Mampe et al. 
[26] claimed to measure a mean life of 877.6 ± 3 s (remarkably consistent with the current 
world average). Subsequently, it appeared that questions regarding the neutron mean life 
had been resolved, as the mean value varied very little between 1990 and 2010, settling at 
885.6 ± 0.8. However, in 2005, there was already a sign that the mean life was about to 
shift to lower values once again. Serebrov et al. [27] reported a very precise measurement 
of 878.5 ±0.8 which was used in BBN calculations in [28]. This was followed by several 
more recent measurements and reanalyses tending to lower values so that the current world 
average is |29j r n = 880.3 ± 1.1. We will explore the impact of the new lifetime on the light 
element abundances. 

On the observational side, the 4 He abundance determination comes most directly from 
measurements of emission lines in highly ionized gas in nearby low-metallicity dwarf galaxies 
(extragalactic HII regions). The helium abundance uncertainties are dominated by system¬ 
atic effects [301 [32] . The model used to determine the 4 He abundance contains eight physical 
parameters, including the 4 He abundance, that is used to predict a set of ten H and He 
emission line ratios, which can be compared with observations [3BH3B] . Unfortunately, there 
are only a dozen or so observations for which the data and/or model is reliable, and even in 
those cases, degeneracies among the parameters often lead to relatively large uncertainties 
for each system as well as a large uncertainty in the regression to zero metallicity. Newly 
calculated 4 He emissivities [3?J and the addition of a new near infra-red line [36] have led to 
lower abundance determinations [3B[EH], bringing the central value of the 4 He abundance 
determination into good agreement with the SBBN prediction. Moreover, Planck measure¬ 
ments of CMB anisotropies are now precise enough to give interesting measures of primordial 
helium via its effect on the anisotropy damping tail. 

New observations and analyses of quasar absorption systems have dramatically improved 
the observational determination of D/H. Using a a handful of systems where accurate de¬ 
terminations can be made, Cooke et al. m not only significantly lowered the uncertainty 
in the mean D/H abundance, but the dispersion present in old data has all but been erased. 
Because of its sensitivity to the baryon density, D/H is a powerful probe of SBBN and 
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now the small uncertainties in both the data and prediction become an excellent test of 
concordance between SBBN and the CMB determination of the baryon density. 

In contrast to the predicted abundances of 4 He and D/H, the 'Li abundance shows a 
definite discrepancy with all observational determinations from halo dwarf stars. To date, 
there is no solution that is either not tuned or requires substantial departures from Standard 
Model physics. Attempts at solutions include modifications of the nuclear rates mmm or 
the inclusion of new resonant interactions [43HI5]> stellar depletion [46]; lithium diffusion in 
the post-recombination universe 033; new (non-standard model) particles decaying around 
the time of BBN [48H5T] : axion cooling [52] : or variations in the fundamental constants 

[SSI El]. 

In this review, we survey the current status of SBBN theory and its compatibility with 
observation. Using an up-to-date nuclear network, we present new Monte-Carlo estimates of 
the theoretical predictions based on a full set of nuclear cross sections and their uncertainties. 
We will highlight those reactions that still carry the greatest uncertainties and how those 
rates affect the light element abundances. We will also highlight the effect of the new 
determination of the neutron mean life on the 4 He abundance and the abundances of the 
other light elements as well. 

Compatibility with observation is demonstrated by the construction of likelihood func¬ 
tions for each of the light elements [TO] by convolving individual theoretical and observational 
likelihood functions. Our BBN calculations are also convolved with data chains provided 
by the 2015 Planck data release |6]. This allows us to construct 2-dimensional (?/, Y p ) and 
3-dimensional (?/, Y Vl N v ) likelihood distributions. Such an analysis is timely and important 
given the recent advances in the 4 He and D/H observational landscape. In fact, despite the 
accuracy of the SBBN prediction for D/H, the tiny uncertainty in observed D/H now leads 
to a likelihood dominated by theory errors. The tight agreement between D/H prediction 
and observation is in sharp contrast to the discrepancy in 'Li/H. We briefly discuss this 
“lithium problem,” and discuss recent nuclear measurements that rule out a nuclear fix to 
this problem, leaving as explanations either astrophysical systematics or new physics. While 
the 6 Li was shown to be an artifact of astrophysical systematics [55, [56], for now, the 'Li 
problem persists. 

As a tool for modern cosmology and astro-particle physics, SBBN is a powerful probe 
for constraining physics beyond the Standard Model 153- Often, these constraints can be 
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parameterized by the effect of new physics on the speed-up of the expansion rate of the 
Universe and subsequently translated into a limit on the number of equivalent neutrino fla¬ 
vors or effective number of relativistic degrees of freedom, N eS . We update these constraints 
and compare them to recent limits on 7V e ff from the microwave background anisotropy and 
large scale structure. We will show that for the first time, D/H provides a more stringent 
constraint on N e g than the 4 He mass fraction. 

The structure of this review is as follows: In the next section, we discuss the relevant 
updates to the nuclear rates used in the calculation of the light element abundances. We 
also discuss the sensitivity to the neutron mean life. In this section, we present our base¬ 
line results for SBBN assuming the Planck value for 77 , the RPP value for r n and N v = 3. 
In section III, we briefly review the values of the observational abundance determinations 
and their uncertainties that we adopt for comparison with the SBBN calculations from 
the previous section. I 11 section IV, we discuss our Monte Carlo methods, and construct 
likelihood functions for each of the light elements, and we extend these methods to discuss 
limits on N u in section V. The lithium problem is summarized in section VI. A discussion 
of these results and the future outlook appears in section VII. 

II. PRELIMINARIES 

A. SBBN 

As defined in the previous section, SBBN refers to BBN in the context of Einstein gravity, 
with a Friedmann-Lemaitre-Robertson-Walker cosmological background. It assumes the 
Standard Model of nuclear and particle physics, or in other words, the standard set of 
nuclear and particle interactions and nuclear and particle content. Implicitly, this means a 
theory with a number N„ = 3 of very light neutrino flavor^] SBBN also makes the well- 
justified assumption that during the epoch of nucleosynthesis, the universe was radiation 
dominated, so that the dominant component of the energy density of the universe can be 

1 We distinguish between the number of neutrino flavors, three in the Standard Model, and N e g, equal 
to 3.046 in the Standard Model, which corresponds to the effective number of neutrinos present in the 
thermal bath due to the higher temperature from e + e _ annihilations before neutrinos are completely 
decoupled. 


6 



expressed as 


P 


7T~ 


7 7 


- ( 2 + - + -N v ) T 


( 1 ) 


taking into account the contributions of photons, electrons and positrons, and neutrino 
flavors appropriate for temperatures T > 1 MeV. At these temperatures, weak interaction 
rates between neutrons and protons maintain equilibrium. 

At lower temperatures, the weak interactions can no longer keep up with the expansion 
of the universe or equivalently, the mean time for an interaction becomes longer than the 
age of the Universe. Thus, the freeze-out condition is set by 


G 2 F T 5 ~ T wk (7» = H(T f ) ~ C4 /2 T 2 , 


( 2 ) 


where T w k represents the relevant weak interaction rates per baryon that scale roughly as 
T 5 , and H is the Hubble parameter 

H 2 = ^G N p (3) 


and scales as T 2 in a radiation dominated universe. Gp and Gn are the Fermi and Newton 
constants respectively. Freeze-out occurs when the weak interaction rate falls below the ex¬ 
pansion rate, Y w k < H. The /3-interactions that control the relative abundances of neutrons 
and protons freeze out at Tf ~ 0.8MeV. At freeze-out, the neutron-to-proton ratio is given 
approximately by the Boltzmann factor, ( n/p)f ~ e -Am/7y ~ 1/5, where Am = m n — m p 
is the neutron-proton mass difference. After freeze-out, free neutron decays drop the ratio 
slightly to (n/p)bbn — 1/7 before nucleosynthesis begins. A useful semi-analytic description 
of freeze-out can be found in [SHI 53] • 

The first link in the nucleosynthetic chain isp + n—)-d + 7 and although the binding 
energy of deuterium is relatively small, Ep = 2.2 MeV, the large number of photons relative 
to nucleons, // _1 ~ 10 9 causes the so-called deuterium bottleneck. BBN is delayed until 
r]^ 1 exp(—Ep/T) ~ 1 when the deuterium destruction rate finally falls below its production 
rate. This occurs when the temperature is approximately T r\j E b / In?/ 1 Q # \ MeV. 

To a good approximation, almost all of the neutrons present when the deuterium bottle¬ 
neck breaks end up in 4 He. It is therefore very easy to estimate the 4 He mass fraction, 


K 


2 jn/p) 

1 + (n/p) 


« 0.25, 


(4) 
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where we have evaluated Y p using (n/p) ~ 1/7. The other light elements are produced in 
significantly smaller abundances, justifying our approximation for the 4 He mass fraction. D 
and 3 He are produced at the level of about 10~ 5 by number, and 'Li at the level of 1CT 10 
by number. A cooling Universe, Coulomb barriers, and the mass gap at A = 8 prevents the 
production of other isotopes in any significant quantity. 

B. Updated Nuclear Rates 

Our BBN results use an updated version of our code m usi, itself a descendant of 
the Wagoner code ED]. For the weak n p interconversion rates, the code calculates 
the 1-D phase space integrals at tree level; this corresponds to the assumption that the 
nucleon remains at rest. The weak n O p interconversion rates are normalized such that 
we recover the adopted mean neutron lifetime at low temperature and density. To this 
we added order-a radiative and bremsstrahlung QED corrections, and included Coulomb 
corrections for reactions with pe~ in the initial or final states [QT]. We neglected additional 
corrections, because the overall contribution of all other effects is relatively insignificant. 
Finite temperature radiative corrections leads to A Y p ~ 0.0004; corrections to electron 
mass leads to A Y p ~ -0.0001; neutrino heating due to e + e _ annihilation leads to A Y p ~ 
0.0002 [61]. Problems in the original code due to the choice of time-steps in the numerical 
integration [62] have been corrected here. Finally, we included the effects of finite nucleon 
mass [63]] by increasing the final 4 He abundance with AY/ = + 0.0012. We formally adopt 
the Particle Data Group’s current recommended value [29]: r n = 880.3 ± 1.1 for the mean 
free neutron decay lifetime, and assume it is normally distributed. 

In addition to the weak n <->• p interconversion rates, BBN relies on well-measured cross 
sections. The latest update to these reaction rates has been evaluated by the NACRE 
collaboration and released as NACRE-II [2I5J. Only charge-induced reactions were considered 
in NACRE-II, and many more reactions are evaluated there than is relevant for BBN. Those 
reactions of relevance are shown in table [T} The compilation presents tables of reaction rates, 
with recommended “adopted”, “low” and “high” values. We assume the rates are distributed 
with a log-normal distribution with: 



& — In \J Flhigh/ Filowi 


(5) 

( 6 ) 




where Ri ow and Rhigh are the recommended “low” and “high” values for the reaction rates, 
respectively [2Jj • These rates agree well with previous evaluations muzm, largely because 
they depend on the same experimental data. Uncertainties in the NACRE-II evaluation are 
similar to, but tend to be larger than previous studies. This may stem from the method 
used to derive accurate descriptions of the data and the fitting potential models used in 
Distorted-Wave Born Approximation (DWBA) calculations. 

TABLE I. Reactions of relevance for BBN from the NACRE-II compilation. 


d(p, 7 ) 3 He 

d(d, 7 ) 4 He 

d(d, n) 3 He 

d(d,p)t 

t(d, n) 4 He 

3 He(d,p) 4 He 

d(a , 7 ) 6 Li 

6 Li(p, 7 )'Be 

6 Li(p,a) 3 He 

'Li(p, a) 4 He 

'Li(p, 7 ) 8 Be a 



a 8 Be is not in our nuclear network, 8 Be is assumed to spontaneously decay into 2 4 He’s. 

The remaining relevant n-induced rates, p(n, ^)d, 3 He(n,p)f and 'Be(n,p)'Li, need to 
be taken from different sources. We adopt the evaluation from [61] for the key reaction 
p(n, 7 )d. The remaining (n,p) reactions we use are rates taken from [H]. We choose log¬ 
normal parameters in such a way to keep the means and variances of the reaction rates 
invariant. 

C. First Results 

As a prelude to the more detailed analysis given below, we first discuss the BBN predic¬ 
tions at a fixed value of 77 . This benchmark can then be compared to the results of other 
codes. 

Here, we fix ?/ 10 = lO 10 ^ = 6.10. This is related to the value of oob determined by Planck 
[ 6 ] based on a combination of temperature and polarization dats^J The result of our BBN 
calculation at 7710 = 6.10 can be found in Table [n] compared to the fit in [3], based on 
the PArthENoPE code [4j- As one can see, the results of the two codes are in excellent 
agreement for all of the light elements. The results can be quickly compared with the 
2 A straight interpretation of the Planck result based on the TT,TE,EE,+lowP anisotropy data would yield 
7710 = 6.09. However this result already includes the He abundance from BBN. Our choice of r] is discussed 
in Section IV.C and Table IV below. 
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observed abundances given in the next section. However a more rigorous treatment of the 
comparison between theory and observation will be given in section IV. 


TABLE II. Comparison of BBN Results 



Vw 

N v 

Y P 

D/H 

3 He/H 

7 Li/H 

6 Li/H 

This Work 

6.10 

3 

0.2470 

2.579 x 10“ 5 

0.9996 x 10" 5 

4.648X 10“ 10 

1.288 x 10“ 14 

Iocco et al. [3] fit 

6.10 

3 

0.2463 

2.578 x 10“ 5 

0.9983 x 10" 5 

4.646 x 10~ 10 

1.290 x 10“ 14 


III. OBSERVATIONS 

Before making a direct comparison of the SBBN results, we first discuss astrophysical 
observations of light elements. Here we focus on 4 He, D, and the Li isotopes, all of which are 
accessible in primitive environments making it possible to extrapolate existing observations 
to their primordial abundances. BBN also produces 3 He in observable amounts, and 3 He is 
detectable via its hyperhne emission line, but this line is only accessible within Milky Way 
gas clouds that are far from pristine [65]. Because of the uncertain post-BBN nucleosynthetic 
history of 3 He, it is not possible to use these high-metallicity environments to infer primordial 
3 He at a level useful for probing BBN [66]. After an overview of the observational status of 
the four remaining isotopes, we turn to the CMB, which includes constraints not only on 
the baryon density but now also on 4 He and N e g. 

A. Helium-4 

4 He has long since been the element of choice for setting constraints on physics beyond the 
Standard Model. The reasoning is simple: as discussed above, the 4 He abundance is almost 
completely controlled by the number of free neutrons at the onset of nucleosynthesis, and 
that number is determined by the freeze-out of the weak n <(-)■ p rates. The resulting mass 
fraction of 4 He is given in Eq. (|d]). As we have seen, the SBBN result for the Y p dependence 
on the baryon density is only logarithmic and therefore 4 He is not a particularly good 
baryometer. Nevertheless, it is quite sensitive to any changes in the freeze-out temperature, 
Tf, through the relation (J2|. However, strong limits on physics beyond the Standard Model 
m requires accurate 4 He abundances from observations. 
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The 4 He abundance is determined by measurements of He (and H) emission lines in 
extragalactic HII regions. Since 4 He is produced in stars along with heavier elements, the 
primordial mass fraction of 4 He, Y p = p{ A ¥Le)/pi > , is determined by a regression of the helium 
abundance versus metallicity However, due to numerous systematic uncertainties, 

obtaining an accuracy better than 1% in the primordial helium abundance is very difficult 
[301 [SUM]. The theoretical model that is used to extract a 4 He abundance contains eight 
physical parameters to predict the fluxes of nine emission line ratios that can be compared 
directly with observations |^j The parameters include, the electron density, optical depth, 
temperature, equivalent widths of underlying absorption for both H and He, a correction 
for reddening, the neutral hydrogen fraction, and of course the 4 He abundance. Using 
theoretical emissivities, the model can be used to predict the fluxes of 6 He lines (relative 
to H /3) as well as 3 H lines (also relative to H 0). The lines are chosen for their ability to 
break degeneracies among the inputs when possible. For a recent discussion see [381 [39j. 

There is a considerable amount of 4 He data available [33j EH]. A Markov Chain Monte 
Carlo analysis of the 8-dimensional parameter space for 93 H 11 regions reported in [31] was 
performed in [32]. By marginalizing over the other 7 parameters, the 4 He abundance (and its 
uncertainty) can be determined. However, for most the data, the y 2 s obtained by comparing 
the theoretically derived fluxes for the 9 emission lines with those observed, were typically 
very large (^> 1) indicating either a problem with the data, a problem with the model, or 
problems with both. Selecting only data for which 6 He lines were available, and a y 2 < 4, 
left only 25 objects for the subsequent analysis. Further cuts of solutions with for example, 
anomalously high neutral H fractions, or excessive corrections due to underlying absorptions 
brought the sample down to 14 objects that yielded Y = 0.2534 ± 0.0083 + (54 ± 102)O/H 
based on a linear regression and Y p = 0.2574 ±0.0036 based on a weighted mean of the same 
data. 

Recently, a new analysis of the theoretical emissivities has been performed [37]. This 
includes improved photoionization cross-sections and a correction of errors found in the 
previous result. The new emissivities are systematically higher, and for some lines the 
increase in the emissivity is 3-5% or higher. As a consequence, one expects lower 4 He 
abundances using the new emissivities. In [38], the same initial data were used with the 
Below, we will use results based on the inclusion of a tenth line (seventh He line) seen in the near infra-red 

mi¬ 


ll 



same quality cuts, now leaving 16 objects in the final sample. Individual objects typically 
showed 5-10 % lower 4 He abundance yielding Y = 0.2465 ± 0.0097 + (96 ± 122)0/H for a 
regression. Once again, one could argue that the lack of true indication of a slope in the 
data over the restricted baseline may justify using the mean rather than the regression. The 
mean was then found to be Y p = 0.2535 ± 0.0036. The large errors in Y p determined from 
the regression were due to a combination of large errors on individual objects, a relatively 
low number of objects with \ 2 < 4, a short baseline in O/H, and a poorly determined slope 
(though the analysis using the new emissivities does show more positive evidence for a slope 
of Y vs O/H). 

More recently, new observations include a near infrared line at A10830 [36]. The im¬ 
portance of this line stems from its dependence on density and temperature that differs 
from other observed He lines. This potentially breaks the degeneracy seen between these 
two parameters that is one of the major culprits for large uncertainties in 4 He abundance 
determinations. There are 16 objects satisfying y 2 < 6 [39] (there are now two degrees for 
freedom rather than one), with all seven He measured (though these are not exactly the 
same 16 objects used in |38]). Indeed it was found [M] that the inclusion of this line, did in 
fact reduce the uncertainty and leads to a better defined regression 

Y p = 0.2449 ± 0.0040 + (78.9 ± 43.3)0/H (7) 

Unlike past analyses, there is now a well-defined slope in the regression, making the mean, 
Y p = 0.2515 ±0.0017 less justifiable as an estimate of primordial 4 He. The benefit of adding 
the IR He line is seen to reduce the uncertainty in Y p by over a factor of two. This is due to 
the better determined abundances of individual objects and a better determined slope. It 
remains the case, however, that most of the available observational data are not well fit by 
the model. Comparing with the value of Y p given in Table [TT] the intercept of the regression 
(J7]) is in good agreement with the results of SBBN. 


B. Deuterium 

Because of its strong dependence on the baryon density, deuterium is an excellent bary- 
ometer. Furthermore, since there are no known astrophysical sources for deuterium pro¬ 
duction [69] and thus all deuterium must be of primordial origin, any observed deuterium 
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provides us with an upper bound on the baryon-to-photon ratio m EQ. However, the 
monotonic decrease in the deuterium abundance over time indicates that the galactic chem¬ 
ical evolution [75J ESJ affects the interpretation of any local measurements of the deuterium 
abundance such as in the local interstellar medium ra, galactic disk ra. or galactic halo 

m- 

The role of D/H in BBN was significantly promoted when measurements of D/H ratios 
in quasar absorption systems at high redshift became available. In a short note in 1976, 
Adams (771 outlined the conditions that would permit the detectability of deuterium in 
such systems. However it was not until 1997, that the first reliable measurements of D/H 
at high redshift became available ra (we do not discuss here the tumultuous period with 
conflicting high and low measures of D/H). Over the next 20+ years, only a handful of 
new observations became available with abundances in the range D/H = (1 — 4) x 10 -5 
[TUGS!. Despite the fact that there was considerable dispersion in the data (unexpected if 
these observations correspond to primordial D/H), the weighted mean of the data gave D/H 
= (3.01 ± 0.21) x 10~ 5 with a sample variance of 0.68. While the data were in reasonably 
good agreement with the SBBN predicted value (discussed in detail below) using the CMB 
determined value for the baryon-to-photon ratio, the dispersion indicated that either the 
quoted error bars were underestimated and larger systematic errors were unaccounted for, 
or if the dispersion was real, in situ destruction of deuterium must have taken place within 
these absorbers. In the latter case, the highest ratio (~ 4 x 10~ 5 ) should be taken as 
the post-BBN value, leaving room for some post-BBN production of D/H that may have 
accompanied the destruction of 'Li - we return to this possibility (or lack thereof) below. 

In 2012, Pettini and Cooke [80] published results from a new observation of an absorber 
at z = 3.05, with D/H = (2.54 ± 0.05) x 10~ 5 corresponding to an uncertainty of about 2% 
that can be compared with typical uncertainties of 10 - 20 % in previous observations. This 
was followed another precision observation and a reanalysis of the 2012 data along with a 
reanalysis of a selection of three other objects from the literature (chosen using a strict set 
of restrictions to be able to argue for the desired accuracy) [40]. The resulting set of five 
absorbers yielded 

= (2.53 + 0.04) x 10" 5 (8) 

with a sample variance of only 0.05. We will use this value in our SBBN analysis below]/] 

4 Note that the most recent measurement described in [81] has a somewhat larger uncertainty, and its 


13 



C. Lithium 


Lithium has by far the smallest observable primordial abundance in SBBN, but as we 
will see provides an important consistency check on the theory - a check that currently is 
not satisfied! In SBBN, mass-7 is made in the form of stable 'Li, but also as radioactive 'Be. 
In its neutral form, 'Be decays via electron capture with a half-life of 53 days. In the early 
universe, however, the decay is delayed until the universe is cool enough that ' Be can finally 
capture an electron at z ~ 30,000 [82], shortly before hydrogen recombination! Thus 'Be 
decays long after the ~ 3 min timescale of BBN, yet after recombination, all mass-7 takes 
the form of 'Li. Consequently, 'Li/H theory predictions sum both mass-7 isotopes. Note 
also that 'Li production dominates at low 77 , while 'Be dominates at high 77 , leading to the 
characteristic “lithium dip” versus baryon density in the Schramm plot (Fig. [Tj) described 
below. 

A wide variety of astrophysical processes have been proposed as lithium nucleosynthesis 
sites operating after BBN. Cosmic-ray interactions with diffuse interstellar (or intergalactic) 
gas produces both 'Li and 6 Li via spallation reaction such as p cr + 16 O ism —> 6, 'Li + ■ • •, and 
fusion 4 He cr + 4 Hei sm —> 6,1 Li + ■ ■ - [S3]. In the supernova V-process,” neutrino spallation 
reactions can also produce 7 Li in the helium shell via u + 4 He —> 3 He followed by 3 He+ 4 He —> 
' Be + 7 as well as the mirror version of these [833 though the importance of this contribution 
to 'Li is limited by associated n B production [H5]- Finally, in somewhat lower mass stars 
undergoing the late, asymptotic giant branch phase of evolution, 3 He burning leads to high 
surface Li abundances, some of that may (or may not!) survive to be ejected in the death 
of the stars [ 86 ]. Thus, despite its low abundance, 'Li is the only element with significant 
production in the Big Bang, stars, and cosmic rays; by contrast, the only conventional site 
of 6 Li production is in cosmic ray interactions HM. 

To disentangle the diverse Li production processes observationally thus requires measure¬ 
ments of lithium abundances as a function of metallicity. As with 4 He, the lowest-metallicity 
data should have negligible Galactic contribution and point to the primordial abundances. 
To date, the only systems for which such a metallicity evolution can be traced are in metal- 
poor (Population II) halo stars in our own Galaxy. As shown by Spite & Spite [88], halo 
main sequence (dwarf/subgiant) stars with temperatures T eS > 6000 K have a constant Li 

inclusion does not affect the weighted mean in Eq. (|8|) 
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abundance, while Li/H decreases markedly for cooler stars. The hotter stars have thin con¬ 
vection zones and so Li is not brought to regions hot enough to destroy it. These hot halo 
stars that seem to preserve their Li are thus of great cosmological interest. Spite & Spite 
[88] found that these stars with [Fe/H] = log[(Fe/H)/(Fe/H) 0 ] < —1.5 have a substantially 
lower Li content than solar metallicity stars, and moreover the Li abundance does not vary 
with metallicity. This “Spite plateau” points to the primordial origin of Li. Furthermore, 
the Li/H abundance at the plateau gives the primordial value if the host stars have not 
destroyed any of their lithium. 

The 1982 Spite plateau discovery was based on 11 stars in the plateau region. Since then, 
the number of stars on the plateau have increased by more than an order of magnitude. 
Increasingly precise observations showed that the scatter in Li abundances is very small for 
halo stars with metallicities down to [Fe/H] ~ —3 [89T - I93] . 

Recently, thanks to huge increases in the numbers of known metal-poor halo stars, Li 
data has ben extended to very low metallicity. Surprisingly, at metallicities below about 
[Fe/H] < —3, the Li/H scatter becomes large, in contrast to the small scatter in the Spite 
plateau found at slightly higher metallicity. In particular, the trend in extremely metal-poor 
stars is that no stars have Li/H above the Spite plateau value, a few are found near the 
plateau, but many lie significantly below [HU US]- This “meltdown” of the Spite plateau 
remains difficult to understand from the point of view of stellar evolution, but in any case 
seems to demand that at least some halo stars have destroyed their Li. Moreover, as we will 
see below, this possibility has important consequences for the primordial lithium problem. 

While the low-metallicity Li/H behavior is not understood, the Spite plateau remains at 
[Fe/H] ~ —3 to —1.5; lacking a clear reason to discard these data, we will use them as a 
measure of primordial Li. Following ref. [9T] we adopt their average of the non-meltdown 
halo stars having [Fe/H] > —3, giving 

(h) = (1.6 ± 0.3) x Hr 1 " . (9) 

The stars in this sample were observed and analyzed in a uniform way, with Li abundances 
having been inferred from the absorption line spectra using sophisticated 3D stellar atmo¬ 
sphere models that do not assume local thermodynamic equilibrium (LTE). 

Most halo star observations measure only elemental Li, because thermal broadening in 
the stellar atmospheres exceeds the isotope separation between 'Li and 6 Li. However, very 
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high signal-to-noise measurements are sensitive to asymmetries in the A6707 lineshape that 
encode this isotope information. There were recent claims of 6 Li detections, with isotope 
ratios as high as 6 Li/'Li < 0.1 |p?j. The implied 6 Li/H abundance lies far above the SBBN 
value, leading to a putative “ 6 Li problem.” 

However, recent analyses with 3D, non-LTE stellar atmosphere models included surface 
convection effects (akin to solar granulation), and showed that these can entirely explain 
the observed line asymmetry [55l f56j. Thus case for 6 Li detection in halo stars and for a 6 Li 
problem, is weakened. Rather, the highest claimed 6 Li/' Li ratio is best interpreted as an 
upper limit. 


D. The Cosmic Microwave Background 

The cosmic microwave background (CMB) provides us with a snapshot of the universe at 
recombination (^ rec — 1100), and encodes a wealth of cosmological information at unprece¬ 
dented precision. In particular, the CMB provides a particularly robust, precise measure of 
the cosmic baryon content, in a manner completely independent of BBN [6, 1TJT71 1 1 0i2j . Re¬ 
cently, the CMB determinations of N e g and Y p have also become quite interesting. In this 
section we will see how the precision of the CMB has changed the role of BBN in cosmology, 
and enhanced the leverage of BBN to probe new physics. 

The physics of the CMB, and its relation to cosmological parameters, is recounted in 
excellent reviews such as Here, we briefly and qualitatively summarize some of the 

key physics of the recombination epoch (z rec ~ 1100) when the CMB was released. 

Tiny primordial density fluctuations are laid down in the very early universe, by inflation 
or some other mechanism. After matter-radiation equality, z eq ~ 3400, the dark matter 
density fluctuations grow and form increasingly deep gravitational potentials. The baryon- 
electron plasma is attracted to the potential wells, and undergoes adiabatic compression as 
it falls in. Prior to recombination, the plasma remains tightly coupled to the CMB photons, 
whose pressure P 7 oc T 4 acting as a restoring force, with sound speed (? s ~ 1/3. This 
interplay of forces leads to acoustic oscillations of the plasma. The oscillations continue 
until (re)combination, when the universe goes from an opaque plasma to a transparent gas 
of neutral H and He. The decoupled photons for the most part travel freely thereafter until 
detection today, recording a snapshot of the recombining universe. 
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Because the density fluctuations are small, perturbations in the cosmic fluids are very 
well described by linear theory, wherein different wavenumbers evolve independently. Modes 
that have just attained their first contraction to a density maximum at recombination occur 
at the sound horizon ~ c s t rec . This lengthscale projects onto the sky to a characteristic 
angular scale ~ 1°; this corresponds to the first and strongest peak in the CMB angular 
power spectrum. Higher harmonics correspond to modes at other density extrema. The 
acoustic oscillations depend on the cosmology as well as the plasma properties, which set 
the angular scales of the harmonics as well as the heights of the features, as well as the 
correspondence between the temperature anisotropies and the polarization. This gives rise 
to the CMB sensitivity to cosmological parameters. 

Precision observations of CMB temperature and polarization fluctuations over a wide 
range of scales now exist, and probe a host of cosmological parameters. Fortunately for BBN, 
one of the most robust of these is the cosmic baryon content, usually quantified in the CMB 
literature via the plasma (baryon + electron) density p b written as the density parameter 
fl b = pb/Pcrit where the critical density p crit = 3Hq/8ttG, with H 0 = 100 h kms _1 Mpc _1 
the present Hubble parameter. This in turn is related to the baryon-to-photon ratio rj via 


V = 


Pcrit 

(m)n° 


n 


B 


( 10 ) 


where = 2((3)T$/n 2 is the present-day photon number density. The mean mass per 


baryon (m) in eq. (10) is roughly the proton mass, but slightly lower due to the binding 
energy of helium. As a result, the detailed conversion depends very mildly on the 4 He 
abundance and we have 



The CMB also encodes the values of Y p and N e g at recombination. The effect of helium 
and thus Y p is to set the number of plasma electrons per baryon. This controls the Thomson 
scattering mean free path (ctt u e ) _1 that sets the scale at which the acoustic peaks are 
damped (exponentially!) by photon diffusion. The effect of radiation and thus of N eS is 
predominantly to increase the cosmic expansion rate via H 2 oc p. This also affects the scale 
of damping onset, when the diffusion length is comparable to the sound horizon |105j . Thus, 
measurements of the CMB damping tail at small angular scales (high £) probe both Y p 
and N e ff. In fact, the CMB determinations of these quantities are strongly anticorrelated-a 
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higher Y p implies a lower diffusion mean free path, which is equivalent to a larger sound 
horizon ~ / da/aH and thus lower N eS . 

Since the first WMAP results in 2003 [99] sharply measured the first acoustic peaks, 
the CMB has determined 77 more precisely than BBN. Moreover, both ground-based and 
Planck data now measure the damping tail with sufficient accuracy to simultaneously probe 
all of (fl b h 2 , Y p , N eS ) P, M, TO? ]. Of course, BBN demands an essentially unique relation¬ 
ship among these quantities. Thus, it is now possible to meaningfully test BBN and thus 
cosmology via CMB measurements alone! 

We adopt CMB determinations of (h2 b /i 2 , Y p , iV e ff) from Planck 2015 data, as described 


in detail below, MV C We note in passing that these CMB constraints rely on the present 
CMB temperature, T°, which is held fixed in the Planck analysis we use. I 11 the future, 
T° should be varied in the CMB data evaluation, in order to provide the most stringent 
constraint on the relative number of baryons in the universe. 


IV. THE LIKELIHOOD ANALYSIS 

As data analysis improved, theoretical studies of BBN moved toward a more rigorous 
approach using Monte-Carlo techniques in likelihood analyses (10811111] . Thus, in order 
to make quantitative statements about the light element predictions and convolutions with 
CMB constraints, we need probability distributions for our BBN predictions, for the light 
element observations and CMB-constrained parameters. We discuss here how we propagate 
nuclear reaction rate uncertainties into the BBN light element abundance predictions, how 
we determine the CMB-parameter likelihood distributions, and how we combine them to 
make stronger constraints. 


A. Monte-Carlo Predictions for the Light Elements 

The dominant source of uncertainty in the BBN light-element predictions stem from 
experimental uncertainties of nuclear reaction rates. We propagate these uncertainties by 
randomly drawing rates according their adopted probability distributions for each BBN 
evaluation. We choose a Monte Carlo of size N = 10000, keeping the error in the mean and 
error in the error at the 1% level. It is important that we use the same random numbers 
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for each set of parameters (rj, N u ). This helps remove any extra noise from the Monte Carlo 
predictions and allows for smooth interpolations between parameter points. 

For each grid point of parameter values we calculate the means and covariances of the light 
element abundance predictions. We add the 1 /a /N errors in quadrature to our evaluated un¬ 
certainties on the light element predictions. We have examined the light element abundance 
distributions, by calculating higher order statistics (skewness and kurtosis), and by his- 
togramming the resultant Monte Carlo points and verified that they are well-approximated 
with log-normal or gaussian distributions. 

In standard BBN, the baryon-to-photon ratio ( rj ) is the only free parameter of the theory. 
Our Monte Carlo error propagation is summarized in Figure [lj which plots the light element 
abundances as a function of the baryon density (upper scale) and p (lower scale). The 
abundance for He is shown as the mass fraction Y, while the abundances of the remaining 
isotopes of D, 3 He, and 7 Li are shown as abundances by number relative to H. The thickness 
of the curves show the ±lcr spread in the predicted abundances. These results assume 
N v — 3 and the current measurement of the neutron lifetime r n = 880.3 ±1.1 s. 


Using a Monte Carlo approach also allows us to extract sensitivities of the light element 
predictions to reaction rates and other parameters. The sensitivities are defined as the 
logarithmic derivatives of the light element abundances with respect to each variation about 
our fiducial model parameters |112j . yielding a simple relation for extrapolating about the 
fiducial model: 


X< 



( 12 ) 


where represents either the helium mass fraction or the abundances of the other light 
elements by number. The p n represent input quantities to the BBN calculations (p, N v ,r^) 
and the gravitational constant Gn as well key nuclear rates which affect the abundance X t . 


p n fi refers to our standard input value. The information contained in Eqs. (13 17) are neatly 
summarized in Table m 
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baryon-to-photon ratio 77 


FIG. 1. Primordial abundances of the light nuclides as a function of cosmic baryon content, as 
predicted by SBBN (“Schramm plot”). These results assume N v = 3 and the current measurement 
of the neutron lifetime r n = 880.3 ± 1.1 s. Curve widths show 1 — a errors. 
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TABLE III. This table contains the sensitivities, a n ’s defined in Eq. 12 for each of the light element 


abundance predictions, varied with respect to key parameters and reaction rates. 


Variant 

Y P 

D/H 

3 He/H 

7 Li/H 

6 Li/H 

?? (6.1xlO -10 ) 

0.039 

-1.598 

-0.585 

2.113 

-1.512 

N v (3.0) 

0.163 

0.395 

0.140 

-0.284 

0.603 

Gn 

0.354 

0.948 

0.335 

-0.727 

1.400 

n-decay 

0.729 

0.409 

0.145 

0.429 

1.372 

p(n,7)d 

0.005 

-0.194 

0.088 

1.339 

-0.189 

3 He(n,p)t 

0.000 

0.023 

-0.170 

-0.267 

0.023 

'Be(n,p) 7 Li 

0.000 

0.000 

0.000 

-0.705 

0.000 

d(p, 7 ) 3 He 

0.000 

-0.312 

0.375 

0.589 

-0.311 

d(d, 7 ) 4 He 

0.000 

0.000 

0.000 

0.000 

0.000 

'Li(p,a) 4 He 

0.000 

0.000 

0.000 

-0.056 

0.000 

d(a, 7) 6 Li 

0.000 

0.000 

0.000 

0.000 

1.000 

t(a, 7) 7 Li 

0.000 

0.000 

0.000 

0.030 

0.000 

3 He(ct, 7 )'Be 

0.000 

0.000 

0.000 

0.963 

0.000 

d(d,n) 3 He 

0.006 

-0.529 

0.213 

0.698 

-0.522 

d(d,p)t 

0.005 

-0.470 

-0.265 

0.065 

-0.462 

t(d,n) 4 He 

0.000 

0.000 

-0.009 

-0.023 

0.000 

3 He(d,p) 4 He 

0.000 

-0.012 

-0.762 

-0.752 

-0.012 
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B. The Neutron Mean Life 


As noted in the introduction, the value of the neutron mean life has had a turbulent 
history. Unfortunately, the predictions of SBBN remain sensitive to this quantity. This 
sensitivity is displayed in the scatter plot of our Monte Carlo error propagation with fixed 
t] = 6.10 x 1CU 10 in Figure [2j The correlation between the neutron mean lifetime and 4 He 
abundance prediction is clear. The correlation is not infinitesimally narrow because other 
reaction rate uncertainties significantly contribute to the total uncertainty in ^e. 



FIG. 2. The sensitivity of the 4 He abundance to the neutron mean life, as shown through a scatter 
plot of our Monte Carlo error propagation. 


C. Planck Likelihood Functions 

For this paper, we will need to consider two sets of Planck Markov Chain data, one for 
standard BBN (SBBN) and one for non-standard BBN (NBBN). Using the Planck Markov 
chain data m, we have constructed the multi-dimensional likelihoods for the following 
extended parameter chains, base_yhe and base_nnu_yhe, for the plikHM_TTTEEE_lowTEB 
dataset. As noted earlier, we do not use the Planck base chain, as it assumes a BBN 
relationship between the helium abundance and the baryon density. 

From these 2 parameter sets we have the following 2- and 3-dimensional likelihoods 
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from the C\ II1. TIpLA—base.yt^^b? ^p) ^ d ^PLA—base_nnu_yhe(^b 5 ^-p: ^u) • The 2 (1 i 11 K'lisioi 1 H1 
base_yhe likelihood is well-represented by a 2D correlated gaussian distribution, with means 
and standard deviations for the baryon density and 4 He mass fraction 


u b = 0.022305 ± 0.000225 


Y p = 0.25003 ± 0.01367 


(18) 

(19) 


and a correlation coefficient r = cov(a j b , Y p ) / \J\ai(ijj b )\ai{Y p ) = +0.7200. 

The two parameter data can be marginalized to yield 1-dimensional likelihood functions 


for rj. The peak and 1-cr spread in rj is given in the first row of Table |IV[ The following rows 
correspond to different determinations of rj. In the second-fourth rows, no CMB data is used. 
That is, we fix rj only from the observed abundances of 4 He, D or both. Notice for example, 
in row 2, the value for rj is low and has a huge uncertainty. This is due to the slightly low 
value for the observational abundance ([7]) and the logarithmic dependence of Y p on rj. We 
see again that BBN+h^, is a poor baryometer. This will be described in more detail in the 
following subsection. Row 5, uses the BBN relation between rj and Y p , but no observational 
input from Y p is used. This is closest to the Planck determination found in [Ej, though here 
Y p was taken to be free and the value of rj in the Table is a result of marginalization over 
Y p . This accounts for the very small difference in the results for rj: rj 10 = 6.09 (Planck)] 


rji o = 6.10 (Table IV). Rows 6-8 add the observational determinations of 4 He, D and the 
combination. As one can see, the inclusion of the observational data does very little to affect 
the determination of rj and thus we use rj io = 6.10 as our fiducial baryon-to-photon ratio. 

The 3-dimensional base nnu ylie likelihood is not well-represented by a simple 3D cor¬ 
related gaussian distribution, but since these distributions are single-peaked we can correct 
for the non-gaussianity via a 3D Hermite expansion about a 3D correlated gaussian base 
distribution. Details of this prescription will be given in the Appendix. 

The calculated mean values and standard deviations for these distributions are: 


co b = 0.022212 ± 0.000242 (20) 

N eS = 2.7542 ±0.3064 (21) 

Y p = 0.26116 + 0.01812 (22) 

These values correspond to the peak of the likelihood distribution using CMB data alone. 
That is, no use is made of the correlation between the baryon density and the helium 
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TABLE IV. Constraints on the baryon-to-photon ratio, using different combinations of observa¬ 
tional constraints. We have marginalized over Y p to create ID r/ likelihood distributions. 


Constraints Used 

T) x 10 10 

CMB-only 

6.108 ±0.060 

BBN+Yp 

4.87^-46 

BBN+D 

6.180 ±0.195 

BBN+Yp+D 

6.172 ±0.195 

CMB+BBN 

6.098 ± 0.042 

CMB+BBN+Yp 

6.098 ± 0.042 

CMB+BBN+D 

6.102 ±0.041 

CMB+BBN+Yp+D 

6.101 ±0.041 


abundance through BBN. For this reason, the helium mass fraction is found to be rather 
high. Our value of Y p — 0.261 ± 0.36(2cr) can be compared with the value given by the 
Planck collaboration [IJ] of Y p = 0.263ld;37. 

In this case, we marginalize to form a 2-d likelihood function to determine both r/ and 
N eS . As in the 1-d case discussed above, we can determine rj and N v using CMB data alone. 
This result is shown in row 1 of Table |V| and does not use any correlation between rj and Y p . 
Note that the value of N v given here differs from that in Eq. (21) since the value in the Table 
comes from a marginalized likelihood function, where as the value in the equation does not. 
Row 2, uses only BBN and the observed abundances of 4 He and D with no direct information 
from the CMB. Rows 3-6 use the combination of the CMB data, together with the BBN 
relation between r] and Y p with and without the observational abundances as denoted. As 
one can see, opening up the parameter space to allow N v to float induces a relatively small 
drop ij (by a fraction of 1 a ) and the peak for N v is below the Standard Model value of 3 
though consistent with that value within 1 a. 

We note that we have been careful to use the appropriate relation between rj and u>b via 
Eq. [TlJ Also, in our NBBN calculations we formally use the number of neutrinos, not the 
effective number of neutrinos, thus demanding the relation: N e g = 1.015333AC. For the 
2D base yhe CMB likelihoods, we include the higher order skewness and kurtosis terms to 
more accurately reproduce the tails of the distributions. 


24 















TABLE V. The marginalized most-likely values and central 68.3% confidence limits on the baryon- 
to-photon ratio and effective number of neutrinos, using different combinations of observational 
constraints. 


Constraints Used 

q \o 

N„ 

CMB-only 

6.08 + 0.07 

o £?'7+0.30 
z -°'-0.27 

BBN+Yp+D 

6.10 + 0.23 

2.85 + 0.28 

CMB+BBN 

6.08 + 0.07 

2.91 + 0.20 

CMB+BBN+lp 

6.07 + 0.06 

2.89 + 0.16 

CMB+BBN+D 

6.07 + 0.07 

2.90 + 0.19 

CMB+BBN+Yp+D 

6.07 + 0.06 

2.88 + 0.16 


D. Results: The Likelihood Functions 


Applying the formalism described above, we derive the likelihood functions for SBBN and 
NBBN that are our central results. Turning first to SBBN, we fix N v — 3 and use the Planck 
determination of q as the sole input to BBN in order to derive CMB+BBN predictions for 
each light element. That is, for each light element species Ah we evaluate the likelihood 


C(Xi) oc j £pla —base_ yhe(^6} {^i}) ^T) 


(23) 


where £bbn(^; {Ah}) comes from our BBN Monte Carlo, and where we use the q—u >b relation 
in eq. (11). In the case of 4 He, we use only the CMB q to determine the Ah = F^bbn 


prediction and compare this to the CMB-only prediction. 

The resulting CMB+BBN abundance likelihoods appear as the dark-shaded (purple, solid 
line) curves in Figure [3j which also shows the observational abundance constraints (§ |III[ ) in 
the light-shaded (yellow, dashed-line) curves. In panel (a), we see that the 4 He BBN+CMB 
likelihood is markedly more narrow than its observational counterpart, but the two are in 
near-perfect agreement. The medium-shaded (cyan, dotted line) curve in this panel is the 
CMB-only Y p prediction, which is the least precise but also completely consistent with the 
other distributions. Panel (b) displays the dramatic consistency between the CMB+BBN 
deuterium prediction and the observed high- 2 : abundance. Moreover, we see that the D/H 
observations are substantially more precise than the theory. Panel (c) shows the primordial 
4 He prediction, for which there is no reliable observational test at present. Finally, panel 
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(d) reveals a sharp discord between the BBN+CMB prediction for 'Li and the observed 
primordial abundance-the two likelihoods are essentially disjoint. 



Y p 10 5 xD/H 



789 10 11 12 1234567 

10 6 x 3 He/H 10 10 x 7 Li/H 


FIG. 3. Light element predictions using the CMB determination of the cosmic baryon density. 
Shown are likelihoods for each of the light nuclides, normalized to show a maximum value of 1. 
The solid-lined, dark-shaded (purple) curves are the BBN+CMB predictions, based on Planck 
inputs as discussed in the text. The dashed-lined, light-shaded (yellow) curves show astronomical 
measurements of the primordial abundances, for all but 3 He where reliable primordial abundance 
measures do not exist. For ^He, the dotted-lined, medium-shaded (cyan) curve shows the CMB 
determination of 4 He. 


Figure [3] represents not only a quantitative assessment of the concordance of BBN, but 
also a test of the standard big bang cosmology. If we limit our attention to each element 


26 















in turn, we are struck by the spectacular agreement between D/H observations at z ~ 3 
and the BBN+CMB predictions combining physics at z 10 10 and z ~ 1000. The consis¬ 
tency among all three Y p determinations is similarly remarkable, and the joint concordance 
between D and 7 He represents a non-trivial success of the hot big bang model. Yet this 
concordance is not complete: the pronounced discrepancy in 'Li measures represents the 
“lithium problem” discussed below (j|V|. This casts a shadow of doubt over SBBN itself, 
pending a firm resolution of the lithium problem, and until then the BBN/CMB concordance 
remains an incomplete success for cosmology. 

Quantitatively, the likelihoods in Fig. [3] are summarized by the predicted abundances 


Y p = 0.24709 ± 0.00025 

(24) 

D/H = (2.58 ±0.13) x 10~ 5 

(25) 

3 He/H = (10.039 ± 0.090) x 10“ 5 

(26) 

7 Li/H = (4.68 ±0.67) x 10~ 10 

(27) 

log 10 ( 6 Li/H) = -13.89 ± 0.20 

(28) 


where the central value give the mean, and the error the lcr variance. The slightly differences 
from the values in Table [III arise due to the Monte Carlo averaging procedure here as opposed 
to evaluating the abundance using central values of all inputs at a single rj. 

We see that the BBN/CMB comparison is enriched now that the CMB has achieved 
an interesting sensitivity to Y p as well as rj. This interplay is further illustrated in Figure 
|4j which shows 2-D likelihood contours in the (r/, Y p ) plane, still for fixed N v = 3. The 
Planck contours show a positive correlation between the CMB-determined baryon density 
and helium abundance. Also plotted is the BBN relation for Y p (rj), which for SBBN is a 
zero-parameter curve that is very tight even including its small width due to nuclear reaction 
rate uncertainties. We see that the curve goes through the heart of the CMB predictions, 
which represents a novel and non-trivial test of SBBN based entirely on CMB data without 
any astrophysical input. This agreement stands as a triumph for SBBN and the hot big 
bang, and illustrates the still-growing power of the CMB as a cosmological probe. 


27 



5.8 6 6.2 6.4 

Baryon-to-photon ratio, 10 10 x ?7 


FIG. 4. The 2D likelihood function contours derived from the Planck Markov Chain Monte Carlo 
base yhe [113] with fixed N v = 3 (points). The correlation between Y p and rj is evident. The 
3 -a BBN prediction for the helium mass fraction is shown with the colored band. We see that 
including the BBN Y p (jj) relation significantly reduces the uncertainty in rj due to the CMB Y p — r) 
correlation. 


Thus far we have used the CMB 7] as an input to BBN; we conclude this section by 
studying the constraints on rj when jointly using BBN theory, light-element abundances, 
and the CMB in various combinations. Figure [5] shows the rj likelihoods that result from 
a set of such combinations. Setting aside at first the CMB, the BBN+A" curves show 
the combination of BBN theory and astrophysical abundance observations, £bbn+ v ( r l) — 
f £bbn(Vi X) £ o bsP0 dX, with X e (Yj,,D/H). The CMB-only curve marginalizes over the 
PlanckY p values £cMB-oni y (??) = / £pLA-base_ y he(y4>, Yp) dY p where we use the 7j — uj^ relation 
in eq. ( [IT| ) . The BBN+CMB curve adds the BBN Y p (rj) relation. Finally, BBN+CMB+D 







FIG. 5. The likelihood distributions of the baryon-to-photon ratio parameter, 77 , given various 
CMB and light-element abundance constraints. 

also includes the observed primordial deuterium. 

We see in Fig. [5] that of the primordial abundance observations, deuterium is the only 
useful “baryometer,” due to its strong dependence on 77 in the Schramm plot (Fig. [Tj). By 
contrast, 4 He alone offers no useful constraint on 77 , tracing back to the weak Y p {if) trend in 
Fig .0 The CMB alone has now surpassed BBN+D in measuring the cosmic baryon content, 
but an even stronger limit comes from BBN+CMB. As seen in Fig. [4[ this tightens the 77 
constraint due to the CMB correlation between Y p and 77 Finally BBN+CMB+D provides 
only negligibly stronger limits. The peaks of the likelihoods correspond to the values in 
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Table IV and the tightest constraints are all consistent with our adopted central value 
rj = 6.10 x 10~ 10 . 


V. THE LITHIUM PROBLEM 


As seen in the panels of Fig. [3] above, the observed primordial lithium abundance differs 
sharply from the BBN+CMB prediction [19}. This discrepancy constitutes the “Lithium 
Problem”, which was foreshadowed before CMB determinations of g, and has persisted over 
the dozen years since the first WMAP data release. For a detailed recent review of the 
lithium problem, see [114] , Here we briefly summarize the current status. 

The most conventional means to resolve the primordial lithium problem invokes large 
lithium depletion in halo stars [H>j. As noted above (I jlllC ), recent observations of the Spite 
plateau “meltdown” at very low metallicity, [Fe/H] < —3, seem to demand that some stars 
have depleted their lithium [92; |95]. Could the other plateau halo stars have also destroyed 
their lithium? Such a scenario cannot be ruled out, but raises other questions that remain 
unanswered: why is the Li/H dispersion so small at metallicities above the “meltdown”? 
And why is there a “lithium desert” with no stars having lithium abundances between the 
plateau and the primordial abundance? 

It is worthwhile to find other sites for Li/H measurements, as clearly halo star lithium 
depletion is theoretically complex and observationally challenging. Unfortunately, the CMB 
itself does not yet provide an observable signature of primordial lithium [96j. However, 
a promising new direction is the observation of interstellar lithium in low-metallicity or 
high-z galaxies ES33- Interstellar measurements in the Small Magellanic Cloud (metallicity 
~ 1/4 solar) find Li/H ISMgMC = (4.8 ± 1.8) x 1CU 10 [115]. This value is consistent with the 
CMB+BBN primordial abundance, but the SMC is far from primordial, with a metallicity 
of about 1/4 solar. Indeed, the SMC interstellar lithium abundance agrees with that of 
Milky Way stars at the same [Fe/H], which are disk (Population I) stars in which Li/H is 
rising from the Spite plateau due to Galactic production. Thus we see consistency between 
lithium abundances at the same metallicity, but measured in very different systems with 
very different systematics. This strongly suggests that stellar lithium depletion has not 
been underestimated, at least down to this metallicity. Moreover, this observation serves as 
a proof-of-concept demonstration that measurements of interstellar lithium in galaxies with 
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lower metallicities would could strongly test stellar depletion and potentially rule out this 
solution to the lithium problem. 

Another means of resolving the lithium problem within the context of the standard cos¬ 
mology and Standard Model microphysics is to alter the BBN theory predictions due to 
revisions in nuclear reaction rates mmm- But as we have seen, all of the reactions 
that are ordinarily the most important for BBN have been well measured at the energies 
of interest. Typically, cross sections are known to ~ 10% or better, and these errors are 
already folded into Fig. [3} A remaining possibility is that a reaction thought to be unim¬ 
portant could contain a resonance heretofore unknown, which could boost its cross section 
enormously, analogously to the celebrated Hoyle 12 C resonance that dominates the 3a —> 12 C 
rate j 1 16 . 

In BBN, the densities and timescales prior to nuclear freezeout are such that only two- 
body reactions are important, and it is possible to systematically study all two-body reac¬ 
tions that enhance the destruction of 'Be. A small number candidates emerge, for which one 
can make definite predictions of the needed resonant state energy and width: 'Be(d, 7 ) 9 B, 
'Be( 3 He, 7 ) 10 C, and 'Be(f, 7 ) 10 B [131 - H5] . However, measurements in 'Be(d, d) 7 Be [ ITT] . 
9 Be( 3 He, t) 9 B [ 118] . and an R -matrix analysis of 9 B [119J all rule out a 9 B resonance. Sim¬ 
ilarly, 10 C data rule out the needed resonance in 10 C 1120 . The upshot is that a “nuclear 
option” to the lithium problem is essentially excluded. 


It is thus a real possibility that the lithium problem may point to new physics at play 
during or after nucleosynthesis. A number of possible solutions have been proposed and 
are discussed in the reviews cited above. Here we simply note that a challenge to all such 
models is that they must reduce 7 Li substantially, yet not perturb the other light elements 
unacceptably. Generally, there is a tradeoff between 'Be destruction and D production 
(usually as a by-product of 4 He disruption) [76), 1121] . Essentially all successful models drive 
D/H to the maximum abundance allowed by observations. However, the new very precise 


D/H measurements (£ III B) dramatically reduce the allowed perturbations and will challenge 
most of the existing new-physics solutions to the lithium problem. It remains to be seen 
whether it is possible to introduce new physically-motivated perturbations that satisfy the 
D/H constraint while still solving or at least substantially reducing the lithium problem. 
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VI. LIMITS ON N eS 


Before concluding, we consider a one-parameter extension of SBBN by allowing the num¬ 
ber of relativistic degrees of freedom to differ from the Standard Model value of N v = 3 and 
iVeff = 3.046. Opening this degree of freedom has an impact on both the CMB and BBN. In 
Fig. [ 6 j the thinner contours show the 2D likelihood distribution in the ( 77 , N u ) plane, using 
Planck data marginalizing over the CMB Y p . We see that the CMB N v values are nearly 
uncorrelated with 77 . The thicker contours include BBN information and are discussed below. 



FIG. 6 . The 2D likelihood function contours derived from the Planck Markov Chain Monte Carlo 
base_nnu_yhe jl 13] . marginalized over the CMB Y p (points). Thin contours are for CMB data 
only, while thick contours use the BBN Y p (r 7 ) relation, assuming no observational constraints on the 
light elements. We see that that whereas in the CMB-only case N v and 77 are almost uncorrelated, 
in the CMB+BBN case a stronger correlation emerges. 
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Turning to the effects of N v on BBN, eqs. ([TJ) - (J3| show that increasing the number of 
neutrino flavors leads to an increased Hubble parameter which in turn leads to an increased 
freeze-out temperature, Tf. Since the neutron-to-proton ratio at freeze-out scales as (n/p) ~ 
e -A m/T f i pjg} ier jy leads to higher {n/p) and thus higher Y p [122], As a consequence, one 
can establish an upper bound to the number of neutrinos ra if in addition one has a lower 
bound on the baryon-to-photon ratio [25] as the helium abundance also scales monotonically 
with r]. The dependence of the helium mass fraction Y on both 77 and N v can be seen in 
Figure [7] where we see the calculated value of Y for N v — 2, 3 and 4 as depicted by the blue, 
green and red curves respectively. In the Figure, one clearly sees not only the monotonic 
growth of Y with 77 , but also the strong sensitivity of V with N v . The importance of a lower 
bound on 77 (or better yet fixing 77 ) is clearly apparent in setting an upper bound on N v . 
Prior to CMB determinations of 77 , the lower bound on 77 could be set using a combination 
of D and 3 He observations enabling a limit of N u < 4 [ 124 ] given the estimated uncertainties 
in Y p at the time. More aggressive estimates of an upper bound on the helium mass fraction 
led to tighter bounds on N v m Una ra. The bounds on N v became more rigorous when 
likelihood techniques were introduced [511571 nMH25HT25| . 

While the dependence of Y p on N v is well documented, we also see from Fig. [TJ there 
is a non-negligible effect on D and 'Li from changes in N u [5j. I 11 particular, while the 
sensitivity of D to N v is not as great as that of Y p , the deuterium abundance is measured 
much more accurately and as a result the constraint on N u is now due to both abundance 
determinations as can be discerned from Table [V] 


By marginalizing over the baryon density, we can form 1-d likelihood functions for N v . 
These are shown in Figure [ 8 } In the left panel, we show the CMB-only result by the blue 
curve. Recall that this uses no BBN correlation between the baryon density and helium 
abundance. While the peak of the likelihood for this case is lowest of the cases considered 
(N v = 2.67) its uncertainty 0.30) makes it consistent with the Standard Model. The 
position of the peak of the likelihood function is given in Table [V] for this case as well as the 
other cases considered in the Figures. In contrast the red curve shows the limit we obtain 
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FIG. 7. The sensitivity of the light element predictions to the number of neutrino species, similar 
to Figure [lj Here, abundances shown by blue, green, and red bands correspond to calculated 
abundances assuming N v = 2,3 and 4 respectively. 


purely from matching the BBN calculations with the observed abundances of helium and 
deuterium. In this case, the fact that the peak of the likelihood function is at N v = 2.85 
can be traced directly to the fact that the central helium abundance is Y p = 0.2449. Given 


the sensitivity of Y p to N v found in Eq. 13, the drop in N v from the Standard Model value 
of 3.0, compensates for a helium abundance below the Standard Model prediction closer 
to 0.247. Nevertheless, the uncertainty again places the Standard Model within 1 a of the 
distribution peak. The remaining cases displayed (in green) correspond to combining the 
CMB data with BBN. There are 4 green curves in the left panel and these have been isolated 
in the right panel for better clarity. As one can see, once one combines the BBN relation 
between helium and the baryon density, the actual abundance determinations have only a 
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secondary effect in determining N v which takes values between 2.88 and 2.91. Using the 
CMB, BBN and the abundances of both D and ^He, yields the tightest constraint on the 
number of neutrino flavors N v = 2.88 ± 0.16, again consistent with the Standard ModeQ 
It is interesting to note, that because of the drop in Y p in the most recent analysis [39j, the 
95% CL upper limit on N u is 3.20. 

l 

0.8 

0.6 

0.4 

0.2 
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1234562 3 4 5 

Number of Neutrinos, N„ Number of Neutrinos, N„ 

FIG. 8. The marginalized distributions for the number of neutrinos, given different combinations 
of observational constraints. The left panel shows the likelihood function the case where only CMB 
data is used (blue), only BBN and abundance data is used (red) and when a combination of BBN 
and CMB data is used (green). The four green curves are shown again in the right panel for better 
clarity. 

It is also possible to marginalize over the number of neutrino flavors and produce a 1-d 
likelihood function for r]i 0 as shown in Figure [9} In the left panel, the broad distribution 
shown in red corresponds to the BBN plus abundance data constraint using no information 
from the CMB. Here the baryon density is primarily determined by the D/H abundance. 
When the CMB is added, the uncertainty in r] drops dramatically (from 0.23 to 0.06 or 
0.07) independent of whether abundance data is used. The 5 green curves are almost in¬ 
distinguishable and are shown in more detail in the right panel. Once again the peak of 
the likelihood distributions are given in Table [Vj The values of r) are slightly lower than 
the Standard Model results discussed above. This is due to the additional freedom in the 
5 Of all the cases considered, the one that can best be compared with the results presented by the Planck 
collaboration [6] is the case CMB+BBN+D. We find N e s = 2.94 ± 0.38(2tj) while they quote A r e ff = 
2.91 ± 0.37. While we obtain similar results to other cases, direct comparison is complicated not only by 
the slight difference in the r] — Y relation due to different BBN codes, but also by the adopted value for 
primordial 4 He. 35 
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likelihood distribution afforded by the additional parameter, N v . 



FIG. 9. The marginalized distributions for the baryon to photon ratio ( 77 ), given different combi¬ 
nations of observational constraints. 


For completeness, we also show 2-d likelihood contours in the 7710 — N v plane in Fig. 10 


The three panels show the effect of the constraints imposed by the helium and deuterium 
abundances. In the first panel, only the helium abundance constraints are applied. The 
thinner open curves are based on BBN alone. They appear open as the helium abundance 
alone is a poor baryometer as has been noted several times already. Without the CMB, the 
helium abundance data can produce an upper limit on N u of about 4 and depends weakly 
on the value of 77 . When the CMB data is applied, we obtain the thicker closed contours. 
The precision determination of 77 from the anisotropy spectrum correspondingly produces a 
very tight limit in N v . Here, we see clearly that the Standard Model value of N v = 3 falls 
well within the 68 % CL contour. 


The next panel of Figure [TO] shows the likelihood contours using the deuterium abundance 
data. Once again, the thin open curves are based on BBN alone. In this case, they appear 
open as the deuterium abundance is less sensitive to N u , though we do note that the contours 
are not vertical and do show some dependence on N v as discussed above. In contrast to ^He, 
for fixed N v , the deuterium abundance is capable of fixing 77 relatively precisely. Of course 
when the CMB data is added, the open contours collapse once again into a series of narrow 
ellipses. 
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FIG. 10. The resulting 2-dimensional likelihood functions for the baryon to photon ratio ( r]) and 
the number of neutrinos {N v ), marginalized over the helium mass faction Y p , assuming different 
combinations of observational constraints on the light elements. 


The last panel of Figure 10 shows the likelihood contours using both the 4 He and D/H 
data. In this case, even without any CMB input, we are able to obtain reasonably strong 
constraints on both r) and N u as seen by the thin and larger ellipses. When the CMB data 
is added we recover the tight constraints which are qualitatively similar to those in the 
previous two panels. 

Finally, above in Figure [6j we show the 2-dimensional likelihood function using either 
CMB only (the thin outer curves which trace the density of models results of the Monte- 
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Carlo), or the combination of CMB and BBN (tighter and thicker curves). In the latter 
case, no abundance data is used. 


VII. DISCUSSION 

Big bang cosmology can be said to have gone full circle. The prediction of the CMB was 
made in the context of the development of BBN and of what became Big Bang Cosmology 
USD]. Now, the CMB is providing the precision necessary to make accurate prediction 
of the light element abundances in SBBN. In the Standard Model with N v = 3, BBN 
makes relatively accurate predictions of the light element abundance as displayed by the 
thickness of the bands in Figure [lj These can be compared directly (or convoluted through 
a likelihood function) to the observational determination of the light clement abundances. 
The agreement between the theoretical predictions and the abundance D/H is stunning. 
Recent developments in the determination of D/H has produced unparalleled accuracy [40]. 
This agreement is seen instantly when comparing the likelihood functions of the observations 
with that of the predictions of BBN using CMB data as seen in the second panel of Figure 
[3| The helium data has also seen considerable progress. New data utilizing a near infrared 
emission line [36] has led to a marked drop in the uncertainty of the extrapolated primordial 
4 He abundance [39]. While the error remains large compared with the precision of the BBN 
prediction, the agreement between theory and observation is still impressive. 

Is two out of three okay? Despite the success of the BBN predictions for 4 He and D/H, 
there remains a problem with 'Li [121 IH4j . The predicted primordial abundance is about 
a factor of three higher than the abundance determined from absorption lines seen in a 
population of low metallicity halo stars. The primordial abundance has since 1981 been 
associated with a narrow plateau [88] of abundance measurements. Recently, the extent of 
this plateau has been called into question as a significant amount of downward dispersion 
is seen at very low metallicity ([Fe/H] < —3) [MIES]- If stellar depletion is the explanation 
of the discrepancy between the plateau value and the BBN prediction, it remains to be 
explained why there are virtually no low metallicity stars with abundances above the plateau 
for all metallicities below [Fe/H] < —1.5. If depletion is not the answer, then perhaps the 
lithium discrepancy points to new physics beyond the Standard Model. 

In this review, we have presented the latest combined analysis of BBN predictions using 
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raw CMB data provided by Planck m mg. He have constructed a series of likelihood func¬ 
tion which include various combinations of the CMB, the BBN relation between the baryon 
density and the helium abundance, and various combinations of 4 He and D/H data. We 
presented detailed fits and sensitivities of the light element abundances to the various input 
parameters as well as the dominant input nuclear rates. This allowed us to make relatively 
precise comparisons between theory and observations in standard BBN. The uncertainty in 
the prediction of 4 He remains dominated by the uncertainty in the neutron mean life. We 
also considered a one-parameter extension of SBBN, allowing the number of relativistic de¬ 
grees of freedom characterized by the number of neutrino flavors to differ from the Standard 
Model value of N v = 3. Despite the additional freedom, strong constraints on r) and N v 
were derived. When all abundance data is used in conjunction with BBN and CMB data, 
we obtain a 95% CL upper limit of N u <3.2. As one of the deepest probes in Big Bang 
Cosmology, BBN continues to thrive. 

Going beyond 2015, we expect further improvements in the data which will better test the 
Standard Model. More high resolution data on 4 He emission lines could yield a further drop 
in the uncertainty in primordial helium. One should recall that there are still only a little over 
a dozen objects which are well described by models of the emission line regions. That said, 
there are less than half a dozen quasar absorption systems which yield high precision D/H 
abundances. Moreover, the nuclear physics uncertainties in D/H now dominate the error 
budget. Thus there is strong motivation for future measurements of the rates most important 
for deuterium: d(p } 7 ) 3 He, as well as d(d, n) 3 He, d(d,n)t, and n(p, 7 ) 3 He [131. fl32 j. We can 
be hopeful that future measurements lead to a reduction in the already small uncertainty 
in primordial D/H; futuristically, there is hope of detecting cosmological 92 cm deuterium 
hyperhne lines that would probe D/H at extremely high redshift [133] . Lastly, we can 
be certain to expect updated results from the CMB data when the Planck collaboration 
produces its final data release. 


APPENDIX 

In order to reproduction the likelihoods, we need order 3 and 4 (skewness and kurtosis) 
terms in the multi-dimensional expansion. 
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In ID, this expansion looks like: 


x [1. + SH 3 (z) + KH 4 (z )], (29) 

(x — /i is the mean value and a is the standard deviation and S and K 



(29) 


where z = (x — u)/a, u is the mean 


are the skewness and kurtosis coefficients needed to adequately describe the distribution via 
a Hermite polynomial expansion. The skewness and kurtosis coefficients are proportional 
to the Markov Chain average of the respective Hermite polynomials (S' oc (H 3 (x)) and 
K oc (H-i(z))). In multiple dimensions, the simple gaussian base distribution is replaced 
with the fully correlated multi-dimensional gaussian: 



(30) 
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